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ABSTRACT 

The luminosity function of galaxies is derived from a cosmological hydrodynamic sim- 
ulation of a A cold dark matter (CDM) universe with the aid of a stellar population 
synthesis model. At z = 0, the resulting B band luminosity function has a flat faint 
end slope of a w —1.15 with the characteristic luminosity and the normalization in a 
fair agreement with observations, while the dark matter halo mass function is steep 
with a slope of a k -2. The colour distribution of galaxies also agrees well with local 
observations. We also discuss the evolution of the luminosity function, and the colour 
distribution of galaxies from z = to 5. A large evolution of the characteristic mass 
in the stellar mass function due to number evolution is compensated by luminosity 
evolution; the characteristic luminosity increases only by 0.8 mag from z = to 2, and 
then declines towards higher redshift, while the B band luminosity density continues 
to increase from z = to 5 (but only slowly at z > 3). 

Key words: galaxies: evolution - galaxies: formation - galaxies: fundamental pa- 
rameters - galaxies: luminosity function, mass function - cosmology: theory 



1 INTRODUCTION 

The CDM model provides us with a basis of our under- 
s tanding of cosmic structure formation and galaxy formation 



(Bmmenthal et al. 1984; Davis et al. 1985). Much of the re- 



cent observational evidence points to a CDM universe dom- 



inated by a cosmological constant A 


(|Efstathiou, Suther- 


land, & Maddox 199C|; Ostriker & Steinhardt 1995J Turnei 


and White 1997; Perlmutter et al. 199£; Ricss et al. 19981 


Balbi et al. 200(4 Lange et al. 2001; 


Hu et al. 2001). The 



physics, of galaxy formation, however , is substantially more 



An alternative approach is to use cosmological hydro- 
dynamic simulations directly. The advantage is that many 
physical processes are automatically taken into account with 
fewer model assumptions. The disadvantage, on the other 
hand, is that the method is computationally expensive, and 
the resolution is limited by computing power. However, the 
Eulerian hydrodynamic mesh is now approaching (1000) 3 , 
and we are perhaps beginning to obtain meaningful results 
on the global properties of galaxies. 

A daptive Mesh Refinement (e.g 



rnmplirlatpH than thp fnrtnatinn af largo-scalp stmctiirp fnr 



Bryan & Norman 



1995; Kravtsov, Klypin, fc Khokhlov 1997) and Smoothed 



which an accurate treatment of gravity suffices , Many 



Parti cle Hydrodynamics (e.g., Katz, Weinberg, fc Hernquist 



thors nbp thp sn-rallpH spmi-analyt.ir mnHpls in which thp 

dark matter (DM) halo formation is supplemented with sim- 
ple models of dissip ative physics for baryo ns. Following the 
pionee ring work of IWhitc fc Frenk (199l|), a series of the 



work (iBaugh et al 1998 



Kauffmann et al. 1998 



somervilk 



and Primack 1998] ; Cole et al. 2000) has shown that the 
model can account for many of the observed galaxy proper- 
ties. 



1996 ) have virtues intermediate between the two listed ap- 
proaches. They solve the hydrodynamic equations directly 
and have a spatial resolution better than the Eulerian 
scheme, but they have a mass resolution coarser than the 
Eulerian scheme adopted in this paper. Semianalytic mod- 
els which are based on dark matter halo merger trees in 
N-body simulations also suffer from mass resolution limits. 



In a precedin g publication (Nagamine, Fukugita, Cen 
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& Ostriker 2001, hereafter Paper I), we discussed galaxy 
formation history, with a special emphasis given to the star 
formation history and the stellar metallicity distribution. In 
this article, we focus our attention to the luminosity func- 
tion (LF) and its evolution, which allows comparison with 
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one of the most fundamental aspects of galaxy observations. 
We supplement our cosmological simulation with a popul a- 
tion synthesis model GISSEL99 ( [Bruzual fc Chariot 199o| ; 
Chariot 1999, private communication), which takes metal- 
licity variations into account. We assume that stars form 
from gas as soo n as the cooling conditions are satisfied. The 
Salpeter (1955) initial mass function (IMF) with a turnover 
at low-mass end reported by Gould, Bahcall, & Flynn (1996 



hereafter GBF) is assumed. The results are presented in the 
standard Johnson-Morgan photometric system. 



2 SIMULATION AND PARAMETERS 

The hydrodynamic cosmological simulation we use in this 
paper is the same as that used in Paper I; the comoving box 
size is 25/i _1 Mpc, with 768 3 grid cells and 384 3 dark matter 
particles each of mass 2.03 x 1O 7 /i _1 M . The comoving cell 
size is 32.6h~ 1 kpc, and the mean baryonic mass per cell is 
3.35 x 10 5 Ii~ 1 Mq. The cosmological parameters are chosen 
to be (n m ,Q A Mbh 2 ,h,a 8 ) = (0.3,0.7,0.016,0.6 7,0.9) . The 



Cen & Os- 



basic s tructu re of the code is similar to that of 

1992a|,|b|), but significantly improved over the years 



trikcr 



(see pen fc Ostriker 2000| ). 

We also use a simulation with 384 3 grid cells (factor 
2 and 8 lower spatial and mass resolution, respectively) to 
assess the effect of the resolution on our result. We call the 
higher resolution run as "N768" , and the other "N384" . The 
input parameters are kept to be same for the two runs except 
the resolution. 

The details of the star formation recipe, the energy feed- 
back, and the metal production treatments are described in 
Paper I, so we do not re peat them here. The adopted yield 
Y — 0.02 (Arnett 1996) dominantly controls the resulting 
amount of stars (f2» = 0.0052 including the 25% ejected gas 
by supernovae) and the metal density, which are consistent 



with the upper limit of the empirical estimate of Fukugita 
Hogan, & Peebles (199S ). We identify gal axies and DM ha- 



los u sing the HOP grouping algorithm (Eisenstein & Hut 



199S ) with the same threshold parameters as in Paper I. 
We refer to the grouped stellar particles as 'galaxies' here- 
after. We confirmed that all galaxies identified by HOP are 
dynamically stable enough to obtain meaningful results for 
the evolution of individual galaxies (see Paper I). 



3 MASS AND LUMINOSITY FUNCTION AT 

Z = 

In Figure [j], we show the DM halo mass function <I>(M) 
in dash-dotted (N768) and short-dashed (N384) histograms. 
The galaxy stellar mass function <E> 9 (M) is also shown in 
solid (N768) and dotted (N384) histograms at z = 0, where 
$(A/)dlog(Af) = 4>(M)dM (likewise for the LF below). The 
two long-dashed lines show 4>(M) ~ M~ 2 and cf> g {M) ~ 

The DM halo distribution follows <j>(M) ~ M~ 2 well, 
which is the generic mass distribution of the hierarchical 
clustering model. Some rounding at the low mass end is 
seen below lO lo A/0 (500 DM particles), and a slight over- 
production of massive halos is seen at above 1O 13 M0 due to 
overmerging problem. 



a 
'2 




Figure 1. Mass function of dark matter halos (dash-dotted N768 
& short-dashed N384 histogram) and galaxy stellar mass (solid 
N768 & dotted N384 histogram). The two long-dashed lines show 
<j}(M) ~ M~ 2 and 4> 9 {M) ~ M -1 - 15 . 



The galaxy stellar mass function behaves quite differ- 
ently: flattening below 10 10 M(Ti is apparent, which can be as- 
cribed to superno va feedback dPekel fc Silk 198(1 ) and pho- 
toheating of gas (Rees 1986; Efstathiou 1992), both more 
efficient in less massive galaxies. The results of N768 and 
N384 agree well in the overlapping mass range for both the 
dark matter and the galaxy stellar mass functions, verify- 
ing that the results that concern us in this paper are not 
strongly affected by the resolution effects. Therefore, only 
N768 run is used in the rest of this paper. 

At above 3 x 10 10 Mq, the galaxy stellar mass func- 
tion departs from the power law and seems to follow DM 
halo mass function. This departure from the faint-end power 
law can partly be ascribed to the inefficient cool ing in mas- 
sive galaxies (|Rees fc Ostriker 1977|; |Silk 1977[). The cut- 



off at high-mass end, however, is not quite exponential; 
there are several objects with a very large stellar mass of 

> 5 x 10 h~ Mq (indicated by the dotted histogram), due 
to an overmerging problem. These objects have the mass of 
groups and clusters, and the galaxies have overmerged in 
very high-density regions due to the limited spatial resolu- 
tion of the simulation. There is no good way to decompose 
the overmerged objects into individual galaxies, and it is dif- 
ficult to estimate the effect of the overmerging to the shape 
of the high-mass end of the stellar mass function. 

The LF in the rest-frame B band (before dust ex- 
tinction) is presented in Figure | as a function of abso- 
lute magnitude (the overmerged objects with stellar mass 

> 5 x 10 11 h~ 1 M@ are indicated by the dotted histograms). 
We attempt to fit the Schechter function to the computed 
LF in the following manner: the faint end slope a is fixed 
to —1.15 throughout this paper, we then choose 4>* to 
fit the faint end plateau of the LF, and adjust Mg so 
that the chosen Schechter function reproduces the com- 
puted luminosity density when it is integrated over. Al- 
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N gal =2097 
z = 0.0 



0* = O.83xlO- 2 h 3 Mpc" 
M*=-21.43 




■25 



Figure 2. Luminosity function of galaxies at z = 0. The galax- 
ies with stellar mass of > 5 X 10 n /i _1 Af o are shown as dotted 
histogram, and the solid curve is the Schechter function with the 
indicated parameters in the figure. See text for the method of the 
fitting. 



3 
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Figure 3. Luminosity function of galaxies at z = 0.5, 1, 3, and 5. 
The 'galaxies' with stellar mass > 5 X 10 11 h —1 Mq are shown as 
dotted histograms. The solid curves are the Schechter functions 
with the indicated parameters in each panel, and the dashed curve 
is the Schechter fit to the simulated LF at z = shown in Figure El 



though the overmerging problem hampers accurate deter- 
mination of the Schechter parameters of the computed LF, 
they can still be used to characterize our result with a min- 
imum number of adjustable parameters. In the figure, the 
Schechter function with a = —1.15, M B = —21.43, and 
4>* = 0.83 x 10~ 2 /i 3 Mpc" 3 is shown in the solid curve. The 
simulated LF fits the slope of a — —1.15 well at fainter mag- 
nitudes below Mb = —17. If the overluminous objects are 
not included in the fit, then we obtain M B — —20.68, which 
is 0.75 magnitude dimmer than the case of full luminosity. 
This may suggest that our M B at z < 1 is biased towards 
the brighter side because of the overmerged objects. 

The computed B band luminosity density Lb of galax- 
ies at z = is 2.4 x 10 8 /iZ/q,b Mpc -3 including the over- 
merged objects. A source of large uncertainty in the lumi- 
nosity density is the assumed IMF. The use of Scalo IMF 
instead of the Salpeter IMF would increase the computed 
Lb by ~ 30% ([White fc Frenk 199l[). Another source of un- 



certainty is the dust extinction. The empirical aver age value 
of extinction in the B band is Ab = 0.33 mag ( [ie Vau" 



couleur ;, et al. 1991 ). The corresponding decrease in Lb by 
dust extinction is by 30%. 

Optical surveys of galaxies have made significant 
prog ress over the last 1 years. In particular, the 2dF Sur- 
vey dCross et al. 200C) and the Sloan D igital Sky Survey 



(Blanton et al. 2001 



Yasuda et al. 2001) have secured the 
luminosity functions and the local luminosity density of the 
universe. Both surveys agree on a = —1.20 ± 0.10 and 



Lb = 2.4 ± 0.4 x 10 8 hL Q>B - Our M B at z = is brighter 
than empirical values by 0.76 magnitude, and <f>* is lower 
by a factor of ~ 2. Overmerged objects contribute about 
50% of the total luminosity in the box, but the overmerging 
does not affect the total luminosity density, as the fraction 



of baryons that condense into stars is basically determined 
by the yield parameter in the simulation (see Paper I) . 

The uncertainty in the fitted value of M B associated 
with the overmerging, IMF, and dust extinction is of the 
same order as the above discrepancy between the computed 
and the observed LF. In view of the fact that we have not 
fine-tuned the simulation parameters, we consider the agree- 
ment of the computed LF and Lb with observations is fair. 



4 EVOLUTION OF LUMINOSITY FUNCTION 

Figure ^ shows the rest-frame B band LF at z — 0.5, 1, 3, 
and 5. The solid curves are the Schechter functions with the 
indicated parameters (see also Table 1), chosen in the same 
manner as in Figure 0. The dashed curve is the Schechter 
function at z = for comparison. 

We observe some evolution between 2 = and 1. The 
characteristic magnitude Mg brightens by 0.8 mag, and the 
normalization <j>* increases by a factor of 1.4. As a result, 
Lb at z = 1 is 2.9 times higher than at z — (Table 1). The 
faint end slope evolves very little. 

We find that the evolution of the LF is different from 
that of the galaxy stellar mass function, which is also fit- 
ted well by the Schechter function except at the high-mass 
end. The characteristic mass M* of the galaxy stellar mass 
function decreases by a factor of 2.4 from z = to 1, but 
the number of galaxies increases by a factor of 1.4, leading 
to a moderate decrease in the stellar mass density p, by a 
factor of 1.5 (Table 1). These changes are compensated by 
the luminosity evolution: the average B band mass-to-light 
ratio of galaxies decreases by a factor of 2.4, leading to more 
luminous L* B at higher redshift. The increase of luminosity 
to higher redshift is due to increased star formation of blue 
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galaxies (Paper I), but it is also partly due to the passive 
evolution of red galaxies, as evidenced from the shift of the 
red edge in the colour distribution, as we see below. 

Empirical knowledge of the evolution of the LF is still 
controversial, despite much eff ort made by th e Canada- 
France Redshift Surv ey (CFRS; [Lilly et al. 1995| ), the Aut- 
ofib Redshift Su rvey flEllis et al. 19961), the CNOC2 S urvey 
( |Lin et al. 199$ ), th e CADIS Survey (Fried et al. 20o"l| ), and 



the DEEP project ( pavis, et al. 200q ). The CFRS and the 
CADIS claim that the blue population brightens by about 
1 magnitude from z = to 0.8, while the characteristic lu- 
minosity changes little for the red population. On the other 
hand, the CNOC2 and the DEEP claim brightening of char- 
acteristic luminosity for red (or early) populations. The Aut- 
ofib indicates little evolution of L* B for total luminosity func- 
tion, while it claims steepening of the faint end slope. 

Overall, we summarize that the change in L* B from z — 
to 1 deduced from the observations for the total LF is 
modest, no more than ~ 1 mag in the A cosmology, whereas 
the appreciable increase in luminosity density is common to 
all surveys. The values of luminosity density are discrepant 
among authors, and their low-redshift values are also not 
consistent with what we referred to for z = above. By 
compiling all results, we conclude that the increase of Lb 
from z = to 0.5 is by a factor of 1.4 — 2, and by a factor 
of 2 — 3 from z = to 1. The result from our simulation is 
consistent with these data at the upper end for both L* B and 
Lb- We have not seen steepening of the total LF towards 
higher redshift in our simulation. 

The trend of evolution is somewhat different when one 
goes to higher redshift. The characteristic mass of the galaxy 
stellar mass function decreases more rapidly from z = 1 
to 3 (by a factor of 6). This change is faster than that of 
the mass-to-light ratio, resulting in a moderate decrease of 
the characteristic luminosity. The normalization of the LF 
increases rapidly towards higher redshift, but the increase 
in Lb is smaller above z = 2, because galaxies become less 
massive and dimmer at the same time. 

We find that the [/-band luminosity density continues 
to increase from z = to 5, whereas the if-band luminosity 



density gradually decreases towards high-redshift, at least 
from z 2 (some increase is observed trom z — U to 1). 
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Figure 4. Rest-frame B — V colour distribution of galaxies at 
various rcdshifts. See text for discussions. 



the distribution on the red side does not change very much 
as a function of redshift. Another conspicuous fact as one 
goes to higher redshift is the increasingly longer tail towards 
the blue side, indicating star formation in larger fractions 
of galaxies. At high-redshift, this star-forming population 
rapidly increases, and at z = 4 the bulk of the population 
is star-forming. In addition to the increasing star formation, 
the stellar population tends to be bluer due to lower metal- 
licity at high-redshift (see Paper I). 

The U — B colour at z = ranges from —0.05 to 0.4 
in good agreement with local observation s, except that it 
is a little narrower t han the observations ( Buta et al. 1994 ; 
Fukugita et al. 1995). The evolution is qualitatively similar 
to B — V . Assuming a typical e xtinction parameter of Ry = 
3.1 and the extinction law of |Cardelli, Clayton, fc Mathis 
(1989p , one finds that the U — B colour becomes redder by 



GA 



LAXY COLOURS AND THEIR 



0.06 mag 

The V — K colour at z = ranges from 1.9 to 2.9. 
This is perhaps bluer by 0.1 mag at the red edge, and 0.3 
mag bluer at the blue edge than the real world (Fioc & 



EVOLUTION 



Figure ^ presents the B — V colour distribution of galaxies 
at various redshifts. At z = 0, the B — V colour ranges from 
0.55 to 0.95 with a tiny fraction of galaxies (0.6% in num- 
ber, but 12% in luminosity) bluer than 0.5 that represents 
very active star-forming galaxies. The median B — V is 0.70. 
These results are in global agreement with observations for 
local galaxies, which ranges from B — V = 0.45 (the bluest 
Im) to 1 (giant E); B — V = 0.7 is the median colour of Sb 



galax ies (Buta et al. 1994; Fukugita, Shimasaku, & Ichikawa 



1995). In finer details, the calculated distribution is slightly 



narrower than the observations in both blue and red edges. 

As one goes to higher redshift, the red edge is shifted, 
which is ascribed to the younger age of galaxies, since the red 
edge is mainly determined by the maximal available time for 
passive evolution of galaxies up to each epoch. The shape of 



Rocca-Volmerange 1999). This bluer colour is ascribed to 
the GISSEL99 model, with which V - K colour 1 Gyr after 
the burst is only 2.0, and that at 14 Gyrs is 3.0 for a galaxy 
with 40% solar metallicity (see Paper I for the metallicity of 
galaxies in our simulation). 



6 DISCUSSION 

Our ACDM simulation reproduces the global statistics con- 
cerning the luminosity and colours of galaxies fairly well at 
z = considering the uncertainties in the model and the 
simulation. Even the current level of agreement is a signif- 
icant improvement over the previous simulations, and we 
are satisfied with the qualitative agreement for the moment, 
given that we have not done the fine-tuning of the input 
parameters. 
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Table 1. GALAXY STATISTICS PARAMETERS 



Quantities 


2 = 


2 = 0.3 


2 = 0.5 


2 = 1 


2 = 2 


2 = 3 


2 = 5 


M B 


-21.43 


-21.93 


-22.21 


-22.24 


-22.27 


-21.94 


-21.12 


<t>* 


0.83 


0.93 


0.97 


1.15 


1.60 


2.50 


5.50 


C B 


1.62 


2.86 


3.86 


4.69 


6.74 


7.77 


8.02 


(M/ L b ) 


4.3 


3.2 


2.7 


1.8 


0.87 


0.44 


0.16 


(M/L K ) 


1.3 


1.0 


0.91 


0.69 


0.42 


0.30 


0.21 


p*/£b 


4.0 


2.1 


1.4 


0.94 


0.44 


0.25 


0.13 


P*/C k 


0.77 


0.53 


0.49 


0.39 


0.27 


0.21 


0.15 


(B-V) 


0.70 


0.68 


0.66 


0.61 


0.50 


0.36 


0.12 


(U-B) 


0.09 


0.07 


0.06 


0.05 


0.03 


-0.03 


-0.28 


(V-K) 


2.29 


2.23 


2.18 


2.09 


1.92 


1.68 


1.20 


m; 


2.34 


1.97 


1.45 


0.97 


0.41 


0.16 


0.05 


*1 


0.83 


0.93 


1.10 


1.35 


2.20 


3.60 


6.00 


p* 


6.50 


6.13 


5.34 


4.39 


3.00 


1.98 


1.03 



Parameters of the Schechter functions shown in Figure g and |3j are presented in the first two rows, a is fixed 
to —1.15 in all cases. See text for the method of choosing the values of Mg [mag ] and </>*[10- 2 h 3 Mpc" 3 ]. 
Cb is in units of [10 8 Lq,b Mpc~ 3 ] with h = 0.67. The average mass-to-light ratio of all galaxies in solar 
units, and the average of the entire box (mass density divided by luminosity density) are also given. The 
Schechter parameters of the galaxy stellar mass function (M*[10 n M 9 ] and <j>* [10" 2 fe 3 Mpc -3 ]) and the 
comoving stellar mass density p»[lO 8 Af0 Mpc — 3 ] are shown in the bottom three rows. 



The number evolution and the luminosity evolution of 
galaxies takes place at the same time in the simulation. The 
galaxies in the simulation continues to become less massive 
and more numerous from z — to 5, as expected in the 
hierarchical structure formation scenario, but luminosity per 
stellar mass significantly increases to high redshifts. Because 
of the compensation of the two effects, the change in the 
LF is modest: the characteristic luminosity increases from 
z = to 2 only by 0.8 mag, and then decreases beyond 
z > 2. The evolution of the LF is a modest representation of 
the evolution of individual galaxies. The B band luminosity 
density continues to increase from z — to 5, but only slowly 
at z > 3. We do not see the apparent steepening of the faint 
end slope of the LF. Uncertainties in the observational LFs 
at non-zero redshifts do not allow a finer comparison of the 
model prediction with the data. 

In detailed level the agreement between the predictions 
and the observations is not perfect. However, given the un- 
certainties in the theoretical modelling and observations, the 
level of agreement presented in this paper is encouraging. 
The qualitative evolutionary trends presented in this paper 
(and those in Paper I) can be taken as predictions of hydro- 
dynamic simulations based on the ACDM model plus the 
prescriptions for star formation and feedback effects. These 
qualitative predictions can be used to test the ACDM sce- 
nario against further observational studies. 



ACKNOWLEDGMENTS 

We thank Michael Strauss for useful comments on the draft. 
K.N. is supported in part by the Physics Department. M.F. 
is supported in part by the Raymond and Beverly Sackler 
Fellowship in Princeton and Grant-in- Aid of the Ministry of 
Education of Japan. R.C. and J. P.O. are partially supported 
by grants AST 98-03137 and ASC 97-40300. 



REFERENCES 

Arnett D., 1996, Supernovae and Nucleosynthesis, Prince- 
ton University Press, pp.496 
Balbi A., et al, 2000, ApJ, 545, LI 

Baugh C. M., Cole S., Frenk C. S., Lacey C. G, 1998 ApJ, 
498, 504 

Blanton M. R., et al., 2001, AJ, 121, 2358 

Blumenthal G. R., Faber S. M., Primack J. R., & Rees M. 

J., 1984, Nat, 311, 517 
Bruzual A. G, Chariot S., 1993, ApJ, 405, 538 
Bryan G. L., Norman, M. L., 1995, AAS, 187.9504 
Buta R., Mitra S., de-Vaucouleurs G, Corwin JR. H. G, 

1994, AJ, 107, 118 
Cardelli J. A., Clayton G. C, Mathis J. S., 1989, ApJ, 345, 

245 

Cen R., Ostriker J. P., 2000, ApJ, 538, 83 
Cen R., Ostriker J. P., 1992a, ApJ, 393, 22 
Cen R., Ostriker J. P., 1992b, ApJ, 399, L113 
Cole S., Lacey C, Baugh C, Frenk C, 2000, MNRAS, 319, 
168 

Cross N., et al., 2000, MNRAS, 324, 825 

Davis M., Newman J., Faber S., Phillips A., 2000, 

Proc. of the ESO/E CF/STSCI worksh op on Deep Fields 

(Springer), in press ( istro-ph/0012189) 
de Vaucouleurs G. 



et al., 1991, Third Reference Cata 
of Bright Galaxies, Springer- Verlag, p. 31 
Efstathiou G, 1992, MNRAS, 256, L43 
Davis M., Efstathiou G, Frenk C. S., White S. D. M., 1985, 

ApJ, 292, 371 
Dekel A., Silk, J., 1986, ApJ, 303, 39 
Efstathiou G, Sutherland W. J., Maddox S. J., 1990, Nat, 
348, 705 

Eisenstein D. J., Hut P., 1998, ApJ, 498, 137 

Ellis R. S., Colless M., Broadhurst T., Heyl J., Glazebrook 

K., 1996, MNRAS, 280, 235 
Fioc M., Rocca-Volmerange B., 1999, A&A, 351, 869 
Fried J. W., et al., 2001, A&A, 367, 788 
Fukugita M., Hogan C. J., Peebles P. J. E., 1999, ApJ, 503, 

518 



© 2001 RAS, MNRAS 000, [|-| 



6 K. Nagamine et al. 



Fukugita M., Shimasaku K., Ichikawa T., 1995, PASP, 107, 
945 

Gould A., Bahcall J. N., Flynn C, 1996, ApJ, 465, 759 
Hu W., Fukugita M., Zaldarriaga M., Tegmark M., 2001, 
ApJ, 549, 669 

Kauffmann G., Colberg J. M., Diaferio A., White S. D. M., 

1999, MNRAS, 303, 188 
Katz N., Weinberg D. H., Hernquist L., 1996, ApJ, 105, 19 
Kravtsov A. V., Klypin A. A., Khokhlov A. M., 1997, ApJ, 

111, 73 



Lange A. E., et al., 2000, Phys.Rev. D63, in press (astro- 
ph/0 |)05004l ) 

Lilly S. J., Tresse L., Hammer F., Crampton D., Le Fevre 

O., 1995, ApJ, 455, 108 
Lin H., et al., 1999, ApJ, 518, 533 

Nagamine K., Fukugita M., Cen R., Ostiker J. P., 2001, 

ApJ, in press (astro-ph/0011472; Paper I) 
Ostriker J. P., Steinhardt P. J., 1995, Nat, 377, 600 
Perlmutter S., et al., 1998, Nat, 391, 51 
Riess A. G., et al., 1998, AJ, 116, 1009 
Rees M. J., 1986, MNRAS, 218, 25 
Rees M. J., Ostriker J. P., 1977, MNRAS, 179, 541 
Salpeter E. E., 1955, ApJ, 121, 161 
Rieke G. H., Lebolsky M. J., ApJ, 288, 618 
Silk J., 1977, ApJ, 211, 638 

Sneden C, Gehrz R. D., Hackwell J. A., York D. G., Snow 

T. P. 1978, ApJ, 223, 168 
Somerville R. S., Primack J. R., 1999, MNRAS, 310, 1087 
Steidel C. C, Adelberger K. L., Giavalisco M., Dickin- 
son M., Pettini M., 1999, ApJ, 519, 1 
Steidel C. C, Adelberger K. L., Dickinson M., Gi- 
avalisco M., Pettini M., Kellogg M., 1998, ApJ, 492, 428 
Turner M. S., White M. 1997, Phys. Rev. D, 56, 4439 

White S. D. M., Frenk C. S., 1991, Ap J 379, 52 

Yasuda N., et al., 2001, AJ, submitted (|astro-ph/0105545j) 



This paper has been typeset from a TfrjX/ DTgX file prepared 
by the author. 



